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Abstract 

Recent tests aimed at the detection of the general relativistic gravit- 
omagnetic Lense-Thirring effect in the gravitational field of the Earth 
with the LAGEOS and LAGEOS II satellites have proved to be af- 
fected, among other things, by the corrupting impact of the secular 
variations J4 and Jq of the second and third even zonal harmonic co- 
efficients of the multipolar expansion of the terrestrial gravitational 
potential. Unfortunately, they are, at present, poorly known so that 
their impact on the performed measurement is ~ 13% over 11 years. 
Also the impact of the static part of J4 and Jg is relevant inducing a 
systematic bias of 4-9% according to different gravity models. More- 
over, the currently available values for them may retain an a priori 
'memory' effect of the Lense-Thirring effect itself. In this paper we 
suggest a novel method to determine J2, J4, Jq along with their sec- 
ular variations J2, J a, J&- Our approach is based on the use of three 
linear combinations of the nodes of LAGEOS, LAGEOS II and Ajisai 
and the perigee of LAGEOS II. The retrieved values for the even zonal 
harmonics of interest are, by construction, independent of each other 
and of the post-Newtonian precessions. The so obtained mini-model 
could subsequently be used in order to enhance the accuracy and the 
reliability of the measurements of the Lense-Thirring effect by means of 
LAGEOS and LAGEOS II satellites. Moreover, it would also allow for 
a clearer investigation of possible seasonal and interannual variations 
of J a and Jq. 

Keywords: Even zonal harmonics, Secular variations of the even zonal 
harmonics, LAGEOS, LAGEOS II and Ajisai satellites, New Earth gravity 
models, Post-Newtonian Lense-Thirring effect 
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1 Introduction 



The long and accurate records of data available from the laser-ranged geode- 
tic satellites have proven to be an excellent and unique tool for investigating 
the long-term variations of the terrestrial gravity field by means of the satel- 
lite laser ranging (SLR) technique. For example, Cheng and Tapley (2004) 
recently analyzed 28 years of data from 1976 to 2003. 

In particular, the secular variations of the even zonal harmonic coeffi- 
cients Ji of the Newtonian multipolar expansion of the Earth's gravitational 
potential have recently received attention (Cheng et al. 1997; Bianco et al. 
1998; Cox and Chao 2002, Dickey et al. 2002, Cox et al. 2003, Cheng and 
Tapley, 2004). This is mainly due to the observed inversion of the rate of 
change of the Earth's quadrupole mass moment coefficient J2 which, since 
1998, began increasing (Cox and Chao 2002). It is not yet clear if such an 
effect is a long-term feature or is short-term in nature. The variations of the 
even zonal harmonics are mainly related to the Earth's lower mantle viscos- 
ity features (Ivins et al. 1993). Recent studies have shown that also ice mass 
losses, and oceanographic and hydrological contributions are important. 

An accurate knowledge of the secular variations of the other low-degree 
even zonal harmonics, with particular emphasis on J4 and Jq, is of the ut- 
most importance for a reliable and consistent evaluation of the total accuracy 
in the tests aimed at the detection of the post-Newtonian gravitomagnetic 
Lense-Thirring effect (Lense and Thirring 1918). It consists of secular pre- 
cessions affecting the longitude of the ascending node O and the argument 
of perigee u of the orbit of a test particle moving in the gravitational field 
of a central spinning body. They are 
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in which G and c are the Newtonian constant of gravitation and the speed 
of light in vacuum, respectively, S is the proper angular momentum of the 
central body, a, e, i are the semimajor axis, the eccentricity and the inclina- 
tion to the body's equator, respectively, of the test particle's orbit. For the 
LAGEOS satellites, whose orbital parameters are in Tabled they amount to 
a few tens of milliarcseconds per year (mas yr _1 ). The major source of bias 
is represented by the corrupting effects induced by the Newtonian part of 
the terrestrial gravitational field. A strategy to overcome this problem was 
put forth for the first time by Ciufolini (1996) who designed a linear combi- 
nation involving the nodes of LAGEOS and LAGEOS II and the perigee of 
LAGEOS II in order to remove the bias due to the first two even zonal har- 
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monies J2 and J4. In the tests performed with such a combination (Ciufolini 
et al. 1998) and the Earth's gravity model EGM96 (Lemoine et al. 1998) a 
total error of 20% was claimed. However, this estimate is largely optimistic 
mainly due to the impact of the non-gravitational perturbations affecting 
the perigee of LAGEOS II and of the uncancelled even zonal harmonics. 
More realistic evaluations point toward ~ 50 — 80% (Ries et al 2003a; Iorio 
2005a). Ries et al. (2003b) suggested for the first time to use only the nodes 
of LAGEOS and LAGEOS II for measuring the Lense-Thirring effect in view 
of the expected improvements in our knowledge of the Earth's gravitational 
field from the GRACE missions. In (Iorio 2003a; Iorio and Morea 2004; 
2005c) the following combination of the residuals of the nodes of LAGEOS 
and LAGEOS II was explicitly proposed 

8ti h + P i5ti hn , (2) 

with pi = 0.546 (see Section |2J). Eq.@ would entirely absorb the grav- 
itomagnetic signature, which is a linear trend with a slope of 48.1 mas 
yr _1 , because the Lense-Thirring effect is purposely switched off in the dy- 
namical force models used for constructing the (O)-(C) residuals 5£1. The 
coefficient p\ allows to cancel out the aliasing effects due to the static and 
time- varying components of the first even zonal harmonic J2. The combi- 
nation of Eq.@ is, instead, affected by the other higher degree even zonal 
harmonics J4, Jq, Jg.. along with their secular variations J4,Jq- While the 
static parts of Ji induce linear precessions, the shift due to the secular vari- 
ations Ji is quadratic. Ciufolini and Pavlis (2004) adopted the combination 
of Eq.© for analyzing almost 11 years of LAGEOS and LAGEOS II data 
with the 2nd generation GRACE-only EIGEN-GRACE02S model (Reigber 
et al. 2005). They claim a total error of 5% at 1-sigma level. Such estimates 
have been criticized by Iorio (2005a; 2005b) who yields a 19-24% 1-sigma 
total error budget estimate. In (Iorio 2005b) it has been shown that the 
systematic error due to the static part of the even zonal harmonics is still 
model-dependent ranging from 4.3% for the GeoForschungsZentrum (GFZ) 
EIGEN-GRACE02S model to 8.7% for the Center for Space Rese arch (CSR) 
GGM02S model (see http:/ /www.csr. utexas.edu/grace/gravity/ ). These fig- 
ures refer to 1-sigma upper bounds obtained by linearly summing the indi- 
vidual mismodelled precessions. It is important to note that the sensitivity 
of GRACE to the low-degree even zonal harmonics has recently been ques- 
tioned (Wahr et al. 2004); this could prevent from obtaining notable future 
improvements of the systematic error in the Lense-Thirring measurement 
with Eq.@ due to the static part of the even zonals. 
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In regard to the secular variations of the even zonal harmonics, at present 
there is a large uncertainty about the magnitude and even the sign of J4 and 
Jq; see (Cox et al. 2003, Table 1). Thus, their impact over the performed 
measurement has been evaluated 13% over 11 years (Iorio 2005b). 

In this paper, we propose a strategy in order to measure the first three 
even zonal harmonics of the geopotential independently of each other and 
of the relativistic effects by means of suitable linear combinations of the 
orbital residuals of the LAGEOS satellites and of Ajisai. This approach 
could also be useful in better analyzing the temporal features of J4 and J 6 
independently of J2. 



2 The linear combination approach and its main 
features 

To determine the even zonal harmonic Je, it is convenient to consider the 
temporal evolution, averaged over many orbital revolutions, of those Keple- 
rian orbital elements which, under such condition, undergo secular preces- 
sions. They are the node the longitude of the perigee u> and the mean 
anomaly M; for a generic satellite A we will denote any of them as \I/ A . 
Among them, the node £1 is by far the most accurately measurable, mainly 
due to its insensitivity to the non-gravitational perturbations. On the con- 
trary, the perigee u is affected by a host of non-gravitational perturbations 
which, in many cases, are difficult to correctly model so that their impact on 
the recovery of some parameter of interest by means of the perigee cannot 
often be reliably assessed and evaluated (Ries et al 2003a; 2003b). The same 
holds also for the mean anomaly M for which the indirect perturbations on 
the mean motion n induced by the disturbing accelerations affecting the 
semimajor axis a are to be considered as well. 

It is well known that the multipolar expansion of the Earth's geopotential 
in spherical harmonics (Kaula 1966) leads to classical secular precessions of 
the node, the perigee and the mean anomaly which are linear in the even 
zonal harmonics J^. Moreover, the Einstein's general theory of relativity 
(GTR) predicts that additional secular precessions affecting the node and/or 
the perigee are also present. They are the gravitoelectric Einstein secular 
rate of the perigee (Einstein 1915) 

3nGM 

WGE = cMl-e 2 )' (3) 
in which M is the central body's mass and n = \f GMj a 3 is the Keplerian 
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mean motion, and the gravitomagnetic Lense-Thirring rates of the node and 
the perigee. In Table [21 we quote the post-Newtonian secular precessions for 
the nodes of LAGEOS, LAGEOS II and Ajisai and the perigee of LAGEOS 
II. 

Let us suppose, for the sake of clarity, that we want to measure two 
particular even zonal harmonics Ji and in a relativity-free fashion. The 
orbital residuals of a given Keplerian element account for any mismodelled 
or unmodelled physical phenomenon affecting that element. In order to 
illustrate the method, let us consider the orbital residuals ^obs of the rates 
for three satellites denoted as A, B and C and assume that they entirely 
account for one of the three features that we wish to measure separately of 
each other, (i.e., for the relativistic effects, or the even zonal harmonic of 
degree t or I ). In regard to the even zonal harmonics, this means that we 
are using a 'truncated' Earth gravity model in the force model routines of 
the orbit determination system; in principle, it would include all the even 
zonal harmonics except for that one in which we are interested in. We can 
write 

**obs = ^ + ^V/+^TR^GTR + A A , 

^obs = Vpe + V^Jf +*g TR /i GT R + A B , (4) 
in which the coefficients are defined as 

■ C^class / c \ 

= -aJT- (5) 

The coefficients tyj have been explicitly worked out for Vl/ = f2 and ^ = uj 
from I = 2 to I = 20 (Iorio 2003b) ; it turns out that they are functions of the 
semimajor axis a, the inclination i and the eccentricity e of the considered 
satellite: ^fj = ^j(o,, e, i; GM, R), where R is the Earth's mean equatorial 
radius. For example, the coefficient of degree I = 2 for the node is 

while for the perigee we have 

3 n ( R\ 2 / o . -i \ ,„\ 

In TableEJwe quote the coefficients ^.i for the nodes of LAGEOS, LAGEOS 
II and Ajisai and the perigee of LAGEOS II from I = 2 to I = 6. 
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The quantities A include all the other classical effects, of gravitational 
and non-gravitational origin, which affect ^ and which have been included 
in the routines of the orbit determination systems like, e.g., GEODYN II 
(Pavlis et al. 1998) or UTOPIA (CSR) with the level of accuracy (or mis- 
modelling, if you prefer) characteristic of their models. If such force models 
were perfect, such other effects would not affect the residuals and we 
would have A = 0. In particular, we can assume they include the mismod- 
elled secular precessions induced by the remaining even zonal harmonics of 
degree other than I or I , possible seasonal and interannual variations and 
the mismodelled non- gravitational perturbations. The quantity /^gtr is a 
solved-for least square parameter which accounts for the post-Newtonian 
effects. It is in classical mechanics and 1 in GTR. We can consider Eq.@ 
as a non-homogeneous algebraic linear system of three equations in the three 
unknowns Jg, J/ and /^gtr- The square matrix of coefficients is represented 



by 



/ 



*1 
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\i/A 

GTR 
GTR 
GTR 



(8) 



By defining 



A = 54> A 



obs 



A 



A 



0) 



it is possible to obtain 
J i = 
J e > = 

HGTR = 



B = 5*1 - A 



C 



5^ c 



obs 



A + h l B + h 2 C 
A + hB + k 2 C 
A + ciB + c 2 C 



\i/ A 

GTR 



+ Cl^i T R + C2^g T R : 



(10) 
(11) 

(12) 

(13) 
(14) 
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where 




(15) 



(16) 



(17) 



(18) 



(19) 



ci = 



C2 = 




(20) 



This approach was proposed for the first time by Ciufolini (1996) in the 
context of the Lense-Thirring experiment with the LAGEOS satellites. He 
derived a combination involving the nodes of the LAGEOS satellites and the 
perigee of LAGEOS II. It is a particular case of Eq.JHJ) and Eqs.JTHI-lGni in 
Which ^ A = ft LAGEOS , ^ B = C^LAGEOS n = ^LAGEOS H ^ = 2 , / = 4, 

^gtr = ^LT) ^gtr = <^lt- Later, other combinations have been proposed 
(Iorio 2003a; Iorio and Morea 2004; Iorio 2005c; Iorio and Doornbos 2005d). 
For example, the coefficient p\ of the combination of Eq.© is 



It is very important to note that Eas.(|22 [) -(|24 [) allow to measure the 
selected even zonal harmonics Ji and J pi independently of each other and of 
the relativistic effects which are a priori assumed to be valid and included 
in the force models. 

3 Three possible combinations for J2, J4 and Jq 

Applying the considerations outlined in Section|2]to the first three even zonal 
harmonics, it is possible to obtain three linear combinations which allow to 
determine them independently of each other and of the relativistic effects. 
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Note that the retrieved values for are instead correlated in, e.g., (Cheng 
et al. 1997): the correlation coefficients are -0.75 for J 2 and J4, 0.76 for J 2 
and Jq and -0.86 for J4 and J4 and Jq. We will use the LAGEOS satellites 
and Ajisai because of their high altitude which allow to reduce the impact 
of the remaining uncancelled even zonal harmonics J§, Jxo, J12, ■■■ which do 
affect the proposed combinations. The direct inclusion of the other geodetic 
satellites, like Starlette (a = 7331 km, i = 49.8 deg, e = 0.0205), which has a 
rather eccentric orbit, and Stella (a = 7193 km, i = 98.6 deg, e = 0) would 
introduce more even zonals due to their semimajor axes and inclinations, 
thus enhancing the systematic error induced by the uncertainties in such 
other even zonal coefficients. 
The combination for J 2 is 

5fl L + ai 5fl L 11 + a 2 Sn Ayi + a 3 5ij L 11 , (22) 

withai = 2.865127386220304, a 2 = -0.1125697464440513, a 3 = -0.03288044638241112. 
The expected signal affecting Ea.(|22jl is 

X i^J 2 + j 2 t\ + [seasonal, interannual effects] + [5Js, 5Jxo, ...} + [nongrav.], 

(23) 

with X = h L 2 + ai h h 2 n + a 2&2 + a3^ L 2 11 = -1046.48249580121 deg day -1 . 
For a\, a-i, 03, X the values in Tableland Table[21have been used. Ea.(|22|l is 
purposely built up in order to cancel out the classical perturbations induced 
by J4,Jq, along with their temporal variations, and the post-Newtonian 
secular precessions. This feature can easily be checked by calculating Eq.(|22|) 
with the coefficients VI/ ± of the classical precessions of degree t = 4, 6 (Table 
EJ) and with the general relativistic precessions (Table EJ): the result is zero. 
For J4 we have 

fl L + bx5(l L 11 + b 2 Sn Ayi + b 3 5ij L 11 , (24) 

with 61 = 0.7547006806958291, b 2 = -0.04496658671744361, 63 = -0.01501043654953821. 
The expected signal for Ea, ()24|) is 

Y i^J^ + j^tj + [seasonal, interannual effects] + [5Js, 5Jxo, ■■■] + [nongrav.], 

(25) 

with Y = tl\ + 6i$7 L 4 11 + b 2 tl l + b 3 u L 4 11 = -86.5180659848 deg day -1 . For 
bxi b 2 , ^3) Y the values in Table Eland Table El have been used. Ea, ((2"Hl is, by 
construction, independent of J 2 , Jq and relativity. 

The even zonal harmonic coefficient Jq can be determined independently 
of J 2 , J4 and relativity by means of 

5n L + Cl 5tl L 11 + c 2 5tt Ayi + c 3 Sdj L 11 , (26) 
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withci = 5.747950720464317, c 2 = -1.073942838094944, c 3 = -0.02646944285900744. 
It is expected that Ea.([26|) is affected by 

2 \J§ + Jet) + [seasonal, interannual effects] + [5Js,SJiq, ...] + [nongrav.], 

(27) 

with Z = n L 6 + ciO L 6 n + c 2 tl l + c 3 dj L 6 11 = -1022.9930383702 deg day -1 . 
Also in this case the values in Table El and Table have been used for 
calculating c%, C2, C3, Z. 

4 The analysis strategy 

The residuals of the nodes of LAGEOS, LAGEOS II and Ajisai and of the 
perigee of LAGEOS II should be built up by adopting, as usual in precise 
orbit determination process, a complete suite of dynamical force models and 
some background reference model for the Earth's gravitational field complete 
to all available degree and order, apart from just the even zonal harmonic 
coefficients that we are interested in. Alternatively, it would be possible 
to include in the reference models default values Jg for the even zonals 
of interest in order to determine corrections A Ji to them. In this way, the 
residuals should entirely-or partly-account for the linear and quadratic shifts 
induced by that even zonal harmonic, which we are interested in, and for 
the all other physical effects, which we are not interested in, according to 
the level of accuracy of the force models included in the orbit determination 
system. For the background reference model of the gravitational field it 
would be better to use some Earth gravity models derived from non-SLR 
datasets like records from CHAMP and GRACE. 

The so built residuals time series should be analyzed over a time span 
rather long to average out the various time-dependent perturbations, in par- 
ticular those induced by the non-gravitational forces acting on the perigee 
of LAGEOS II. In addition to the expected linear and quadratic integrated 
shifts of Eq. (|23|) . Eq. (|25|) and Eq. (|27|) also seasonal and interannual varia- 
tions could affect the investigated signals. In (Cheng and Tap ley 2004) it 
has been shown that such additional features mainly affect J 2 : our strategy 
could allow to better investigate their effects on J4 and Jq which should be 
less relevant. 

With a quadratic fit QF = clq + bot + c$t 2 it is possible to retrieve both 
Ji (or AJg) and Ji for the considered even zonal harmonic from bo and cq. 
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4.1 Possible sources of errors in J 2 , J4, Je 

The accuracy of the orbit determination process sets the level of the ob- 
tainable observational errors in the even zonal harmonics jj obs ' 1 . The rms 
orbital accuracy deltar for the LAGEOS satellites is < 1 cm; for Ajisai, 
which is more sensitive to the non-gravitational perturbations, we conserva- 
tively assume 10 cm. By assuming 50 ~ 5r/a and 5uj ~ 5r/ea, we obtain 
54 obs) = 2 x 10- 14 ,54 obs) = 9 x 10- 14 ,^4 obs) = 2 x 10- 12 . over a time 
span of, say, one year. 

A major source of systematic error is represented by the impact of the 
mismodelled even zonal harmonics of higher degree which are not cancelled 
out by Eas.(|22 [) -(|26 [) . For example, Ea.()24[) is independent of J2, Jq and rel- 
ativity but is affected by J§, Jiq, J12, ... which are assumed to be modelled 
according to some Earth gravity model solution. Such bias can be evalu- 
ated, in a conservative way, by using the available sigma 5Jg of the variance 
matrix of the adopted background reference model and by summing up the 
individual secular precessions induced by the uncancelled mismodelled even 
zonal coefficients. For example, the systematic error in the measurement of 
J4 with Eq.pil) can be evaluated as 



r (zon) _ ^ e=i 



o L f + hoh 11 + b 2 df + b 3 dj L e 11 



W Y 



SJe 

— (28) 



The impact of the uncancelled even zonal harmonics on the measurement of 
J4 with Ea.(|24|) can analogously be evaluated as 



oh + bxtlh 11 + b 2 of + b 3 co L f 11 



SJe 

— (29) 



v^20 
i(zon) _ Q Z ^= 8 

where T b s is the adopted observational time span in years. The factor 2 in 
Ea. (|29|) comes from the fact that we have compared the integrated linear 

shift due to Jg, J10, ... to the quadratic shift Y {^t 2 ^ over T Q ^ S . 

According to the combined CHAMP+GRACE+terrestrial gravimetry/altimetry 
EIGEN-CG01C Earth gravity model (Reigber et al. 2004), the systematic 
errors due to the uncancelled even zonal harmonics would be 8J^ = 
6 x 10- 13 ,5j| zon) = 3 x 10- 12 ,<5J^° n) = 8 x 10~ 12 . For the secular varia- 
tions of the even zonal harmonics of interest we have^ zon) = Ixl0- 12 /T obs 

yr"\ 5Jt n) = 3 x 10- 12 /T obs yr" 1 , SJ^ = 1 x lO" 11 /!^ yr" 1 . They 
are 1-sigma upper bounds obtained, in a conservative way, by summing up 
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the absolute values of the individual mismodelled precessions according to 
Eq.(j2Hl and Eq.jHJ). 

It is interesting to note that the proposed determination of J4 and Jq, 
which are the most relevant with respect to the Lense-Thirring effect deter- 
mination, would not be biased by a number of effects which, on the contrary, 
should be considered for J2. Among them, there is the semisecular harmonic 
perturbation induced by the 18.6-year tide because its major power is con- 
centrated just in the I = 2 m = constituent. Consequently, it would affect, 
in principle, the recovery of J2. According to recent estimates of the am- 
plitudes of the perturbations induced by it on the nodes of the LAGEOS 
satellites and the perigee of LAGEOS II and of the level of misomodelling in 
the £?2 Love number (1.5% Iorio 2001) in terms of which they are expressed, 
the 18.6-year tide should have an impact on the measurement of J2 of the 
order of 10~ n over one year. However, it might be a rather pessimistic esti- 
mate. Indeed, there are sufficiently long data records to empirically estimate 
the 18.6-year tide (Eanes and Bettadpur 1996) whose action could thus be 
removed from the times series if Ea, (j22|) . Indeed, the available data records 
of LAGEOS and Starlette are nearly 30 years long. 

In regard to other sources of systematic errors, it is important to notice 
that the impact of the non-gravitational perturbations affecting especially 
the perigee of LAGEOS II should be reduced by the fact that the coefficients 
a 3)^3> c 3 with which it enters the three combinations are all of the order of 
10 -2 . In order to get an-order-of-magnitude estimate, some authors claim 
that the non-gravitational perturbations affecting the perigee of LAGEOS 
II could have an impact on the performed Lense-Thirring tests of the order 
of 100% (Ries et al. 2003a; 2003b), i.e. almost 60 mas yr" 1 (Ciufolini 
1996). If so, the bias induced on our proposed estimation of J2, J4 and 
Jq would be of the order of 10 -13 over a time span of, say, one year. In 
regard to Ajisai, its node is not affected by low-frequency perturbations 
which could bias the recovery of the quantities of interest over not too long 
time spans. For example, the most insidious tidal perturbation, i.e. the 
I = 2, m = 1 K\ tide, has the same period of the node which is 0.32 years 
for Ajisai. It is certainly more sensitive than the LAGEOS satellites to 
the non-gravitational perturbations (Sengoku et al. 1995; 1996) mainly 
due to its larger area-to-mass A/M ratio which is 2.7 x 10 -2 m 2 kg -1 : 
for LAGEOS A/M = 6.9 x 10~ 4 . However, their nominal impact can be 
considered negligible, as can be inferred from (Iorio and Doornbos 2005) in 
which it is shown that their influence on a proposed direct measurement of 
the Lense-Thirring effect is small. Moreover, note that the coefficient 62 of 
the combination for J4, which is particularly important for our purposes, is 
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of the order of 10 . 

5 Conclusions 

Motivated by the need of improving the reliability and the accuracy of the 
Lense-Thirring test performed with the node-node LAGEOS-LAGEOS II 
combination of Eq. (j2J) , we have outlined a strategy to determine J% , J4 , Jq 
along with their secular variations J2, J4, Jq, independently of each other 
and of relativity itself. Indeed, one of the major sources of systematic er- 
ror in such a relativistic test is represented by the aliasing classical linear 
and quadratic shifts due to the mismodelling in J4, Jq, J4, J$. According 
to the most recent CHAMP/GRACE Earth gravity models, the systematic 
error induced by the static part of the zonals ranges from 4% to ~ 9%, 
while the impact of their secular variations is 13% over 11 years. The three 
linear combinations of Eas.(|22 [) -(|26 [) built up with the nodes of LAGEOS, 
LAGEOS II and Ajisai and the perigee of LAGEOS II allow to disentangle, 
by construction, the relativistic precessions from those induced by the even 
zonals whose measured values could thus be safely used for reanalyzing the 
Lense-Thirring combination of Eq.(J2J). In regard to J4 and Jq, it should be 
noted that the proposed technique allows for a determination independent, 
by construction, of the static and time- varying parts of J2. Moreover, a re- 
analysis of Eq. (J2J with the values of J4 and Jq determined with the proposed 
approach would not be driven towards the expected result by possible 'im- 
printing' of relativistic effects. The so obtained low-degree minimodel would 
also have an intrinsic value itself, especially in view of the recently discussed 
ability of GRACE to accurately measure the low-degree even zonal harmonic 
coefficients (Wahr et al. 2004). Moreover, a clearer investigation of possible 
seasonal and interannual variations in J4 and Jq could be performed. 
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Tables 

Tabled Orbital parameters of LAGEOS, LAGEOS II and Ajisai. 
Table |2j Post-Newtonian gravitomagnetic (LT) and gravitoelectric (GE) 
secular precessions, in rad s _1 , for the nodes of LAGEOS, LAGEOS II and 
Ajisai and the perigee of LAGEOS II. The perigees of LAGEOS and Ajisai 
are not good observables because of the too small eccentricities of their 
orbits. 

Table El Newtonian secular precession coefficients 4/ ' £ = 9\l/ciass/<9«^; in rad 
s" 1 , for the nodes of (L), (L II) and (A) and the perigee of LAGEOS II 
(£ = 2, 4, 6). The explicit expressions of Ctj and tbj up to degree £ = 20 can 
be found in Iorio (2003). 
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Table 1: 



Orbital element 


LAGEOS 


LAGEOS II 


Ajisai 


semimajor axis a (km) 


12270 


12163 


7870 


eccentricity e 


0.0045 


0.014 


0.001 


inclination i (cleg) 


110 


52.65 


50 



Table 2: 



PN precession 


LAGEOS 


LAGEOS II 


Ajisai 




4.7 x 10~ 15 


4.8 x 10" 15 


1.79 x 10~ 14 






-8.8 x 10~ 15 








5.143 x 10~ 13 





Table 3: 

^j, 1 = 2 £ = A £ = 6 

n L e 6.4393531015 xl0~ 5 2.3720267441 xl0~ 5 4.994251263 xl0~ 6 

n l e 11 -1.1781974640xl0" 4 -8.582111371xl0- 6 7.668773045 xl0~ 6 

-5.72645349511 xl0~ 4 -2.5332565232 x 10~ 5 2.379851554x 10~ 4 

w L f n 8.1596225662xl0- 5 6.0312610809x 10~ 5 5.3628932 x lO" 6 
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